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1 . Introduction 


Detailed Monte Carlo simulations of extensive air showers to be 
detected by the Homestake Surface-Underground Telescope (1) and other 
similar detectors located at sea level and mountain altitudes have 
been performed for 10^-10*° eV primary energies. The results of 
these Monte Carlo calculations will provide an opportunity to compare 
the experimental data with different models for the composition and 
spectra of primaries and for the development of air showers (2), 

In the present paper we report on the results obtained for exten- 
sive air showers generated by 10^-10^ eV primary protons. 


2, Monte Carlo Simulation of the Hadronic Cascade 


The interaction model used assumed inelastic cross sections in- 
creasing with energy as £n*^s, which corresponds to a £n^s rise of the 
inelastic p-p cross section. The x distribution obeys radial scaling 
and is realized through the splitting technique suggested by A,M, 
Hillas, The K/tt ratio is 0.09 and energy independent. 

Transverse momentum distribution is of the form p x exp (-Kp^ ) 


and includes the "sea gull" effect. For x < 0.2, K= 


x + 1/4* 


while 


for x > 0.2 we used a constant exponent of -4.44 for pions. 

The secondary particles are followed until they reach certain 
threshold energy levels: 1.7 TeV for muons, 0.01 and 0.1 TeV for 
electrons, gammas, and 0.1 TeV for hadrons. 


3. Simulation of the Electromagnetic Cascades 

The electromagnetic cascades generated by gammas and electrons 
are propagated to the height of detection by applying the parametrized 
analytical formulas for the longitudinal development and lateral 
distribution of the electron component obtained by Fenyves and Yunn in 
previous Monte Carlo simulations of electromagnetic cascades (3). The 
threshold energy levels for the electromagnetic cascade generating 
gammas and electrons were: 0.01 TeV for 10^ and 10* eV primary 
protons, and" 0*1 TeV for 10^ eV protons. The average number of 
gammas and electrons above the threshold per shower is given in Table 
1. (Muons above 1.7 TeV are also included in the Table.) 
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From the extrapolation of 
the gamma spectra below the 
threshold levels correction 
factors were calculated for 
the missing low energy elec- 
tromagnetic cascades. The correction factors for the total number of 
electrons of the air shower are dependent on the depth and the primary 
proton energy varying between 1.010 and 1.021. 

The geomagnetic effect was included in the calculation by stret- 
ching the east-west axis of each vertical electromagnetic cascade by 
[1 + 0.05(cos X/P)^]^ where X is the geomagnetic latitude and P the 
pressure in atmospheres (4). The individual cascades were then folded 
together to form the electron component of the extensive air shower. 

4 . Results 


Table 1 

E p (eV) Ny N e N y 

1 0 14 361 ±8 3.9 ± 0.3 0.30 ±0.06 

lO 1 * 2612 ±61 31.9 ± 1.1 1.54 ±0.17 

10 16 3015 ±78 34.9 ± 1.2 8.20 ±0.50 


We have run a total of 220 extensive air showers generated by 
10*, 10 15 and 10*° eV primary protons in a standard atmosphere. As 
expected the longitudinal development of these showers could not be 
approximated well with the standard formula used for the electron 
component of electromagnetic cascades generated by single gammas (3) 


N e (E Q ,E,t)= 


A(E,s)0.31 exp [t(l-l .5 £n s)] 


( 1 ) 


where t is the depth measured in radiation lengths (37.1 g/cm^), y = 


£n E J e Q 


( e Q = 81 MeV), 


s = 


3t 


t + 2y (2) 

and A(E,s) is the fraction of electrons having energies larger than 
the electron threshold energy, E. (E=5MeV and A(5MeV ,s)=0.67 were used 
in the present study (3).) 

We could, however, fit N e raax and t max calculated from Eq.(l) 
simultaneously to the corresponding average values obtained for the 
10*, 10 15 and lO 1 ^ eV showers by varying E Q , the energy of the 

cascade generating gamma. The longitudinal development of the simu- 
lated air showers was, however, increasing faster before t max , and 
decreasing slower after t max than the values calculated from Eq.(l). 
This is expected because the electromagnetic component of the air 
shower starts with a large number of gammas, and dies out slower after 
the maximum due to new electromagnetic cascades generated by the 
hadronic core. 

According to this we modified Eqs.(l) and (2) simply by replacing 
t in both equations by 


4 t?;= t + 3 (t max -OtT (3) 

where 3 and y are constants depending on the primary proton energy, 
and the second term increases or decreases the effective t f values as 
compared to t before or after the maximum, respectively. 

The optimum E Q , 3 and Y values obtained for this 

parametrization for E p =10 14 , 10 15 and 10 16 eV are given in Table 2. 
Fig. 1 shows the good agreement between the Monte Carlo simulated 
average values (dots with error bars) and the smooth curves calculated 
by using Eqs.(l) and (2) modified by Eq.(3) with the parameter 
values given in Table 2. 
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Ep(eV) 

E 0 (eV) 

8 

Y 

10* 4 

5.2 . 10 1 ? 

0.0100 

1.00 

10 15 

6.4 . 10 14 

0.0173 

0.75 

10 16 

6.9 . 10 15 

0.0300 

0.50 


Table 2 

Eqs.(l), (2) and (3) can also 
be used to describe moderate- 
ly inclined showers ( 6 < 30° ) . 
Then t in Eq.(3) is the slant 
depth measured from the top of 

* ” the atmosphere. The energy of 

the single gamma, E Q , simulating the longitudinal development of the 
air shower is increasing with E , the energy of the shower as expected 
because of the increasing probability of charged pions and kaons to 
produce nuclear interactions rather than decay. The 3 and Y valufis can 
be easily interpolated for intermediate Ep values when plotted on a 


log-log or lin-log scale vs. E , respectively. ^ 

The lateral distribution of electrons of air showers generated 
14 inl5 an( j eV protons was compared with the modified NKG 


10 1H , 10 

formula 


f(~, ) = C(s) <£, ) S 
r M r M 


2 £ 


+ 1 ) 


s - 4.5 


M 


i (4) 


where r'w is about half of the Moli^re length length,- rj,j, as 
suggested by A. M. Hillas (5), and s is calculated by Eq.(2> modified 
by Eq.(3) with the parameters given in Table 2. A good agreement was 
found for not too large distances from the core, from r“lm to about 
300m, as illustrated in Figs. 2, 3and4 for 850 g/cm 2 , the height of 
the Homestake Telescope where r'^=45m (3). Similarly good agreement 

was found for other heights too. At larger distances from the core, 
however, even the modified the NKG formula fails to^gree with the 
Monte Carlo simulated data. Corrections for this effect will be 
discussed in a forthcoming paper. 

The geomagnetic effect is relatively small at sea level but 
increasing with altitude. The east-west axis of vertical shpwers is 
stretched by a factor of 1.013 at sea level, and 1.019 at 850 g/ cm • 
The relative difference between the f(r/rjj.) values averaged in a cone 
±45° around the east-west axis and in a cone ±45° around the north- 
south axis is increasing with r from about -0.015 at r*lm, to about 
-0.075 at r= 1000m. 
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